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1. Introduction 

Due to their good position resolution, low material budget and low cost, multi-wire proportional 
chambers (MWPCs) are instrumental in high-energy physics. They are currently employed in 
experiments needing high resolution tracking (ALICE-TPC muon spectrometers (LHCb ^) 
and photo-sensitive detectors (HADES-RICH or ALICE-TRD |^), to mention some examples. 
In particular, a new generation of fast MWPC counters aimed at TR detection has been envisaged 
for the CBM experiment at FAIR [^. The system's conceptual design (presented in [^) is based 
on a configuration aiming at a pion suppression better than 100 for momenta higher than p=l.5 
GeV/c with position resolution of the order of 200-300 /xm at incident fluxes up to = 100 
kHz/cm^. These requirements are driven by the identification of and *F' in the di-electron 
decay channel. Populations of and (in particular the ratio of them) are considered as the 
most promising signatures of Quark Gluon Plasma and its evolution [|8|, ^. However, the very 
small production yields (specially in the case of demand very high interaction rates (up to 10^ 
collisions per second), posing challenges for the detector technologies 

In this work the rate capability of classical MWPCs with a gap of h = 3 mm and s = 3,4 mm 
anode wire pitch has been explored, for binary gas mixtures Xe-C02, Ar-C02 and Ne-C02. CO2 
has been chosen as the UV-quencher due to its low chemical reactivity, non-flammability and very 
good ageing properties while we focused on Xe due to its high cross-section for X-ray absorption. 
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The data is interpreted within the theoretical framework developed by Mathieson and Smith ([jl0|], 
[0). 

This document is structured as follows: the setup is described in section ^ once the character- 
istic gain vs voltage curve has been measured (section the behavior of the gain as a function of 
the rate (section ^ can be described within the Mathieson model (subsection ) after beam-size 
corrections are accounted for (subsection 4.2), allowing for a determination of the mobilities of 
the drifting ions (subsection 4.3). An extrapolation to the most common case of uniform illumina- 
tion by minimum ionizing particles (mips), together with a discussion of the results is presented in 
subsection 4.4 and section B|, respectively. 



2. Experimental setup and description of the measurements 

The experimental setup is shown in Fig. |l|. The chamber wires are made of gold-coated tungsten 
of 10 lJ.m radius. The MWPC has a very thin entrance foil that acts as one of the cathode planes 
(25 /xm aluminized kapton). The window dimensions are 6 cm x 8 cm but a slightly smaller 
fiducial area of 5 cm x 5 cm was determined as the region where the gain uniformity was below 
rmsm < 5%. The chamber was irradiated with an X-ray tube, coUimated down to a nominal area 
of approximately A ~0.5 cm^, centered with respect to the afore-mentioned fiducial area. The 
estimation of A was done via simultaneous exposure of 10 Polaroid films attached to the window, 
producing images like the ones shown in Fig. ^ after developing the film. The duration of the 
exposure was chosen as much as possible such that the first films of the stack were over-exposed 
and the last under-exposed. By looking at the beam-profile as a function of the film number the 
beam-quality can be then assesed. The small presence of tails in some cases (Fig. ^, compare 2-up 
and 5-up) was taken into account as a 10% uncertainty in the area determined by this procedure. 

Unless stated otherwise, all the measurements here presented have been performed with an 
X-ray tube operated at a voltage V = 9 kV, placed approximately 10 cm away from the chamber 
By using different filters (0.25-0.5 mm thick Al foils) most of the low-energy bremsstrahlung and 
characteristic radiation can be suppressed. The feature-less bremsstrahlung spectra can, by means 
of this procedure, render a narrow structure close to the maximum photon energy. 

The measured charge was calibrated by using as reference the internal conversion line of a 
Fe^^ source (Ex-ray = 5.9 kV). Typical charge spectra measured with the ADC after calibration are 
shown in Fig. ^ for a Ar-CO2(80-20) mixture, obtained with the source (squares) and with the tube 
(triangles, circles) for two different filters. The energy spectra obtained with a Fe^^ source showed 
a peak-resolution Ge/E = 10 — 20% with a clear Argon escape peak placed at approximately one 
half of the total absorption one. The spectrum measured with the X-ray tube was broader (Ge/E = 
15 — 30%), resulting from the wider distribution of the bremsstrahlung photons. 

In order to obtain the distribution of Fig. ^ the signal induced in both cathode planes was 
added and amplified by a custom-made charge sensitive pre-amplifier followed by an ORTEC 452 
spectroscopic amplifier. The analog signal was duplicated with a linear Fan-in-Fan-out and sent 
to a LeCroy 0628 leading edge comparator and to a Lecroy 2249W charge ADC. The output of 
the former was used for gating the ADC and for the rate determination, via a scaler The chamber 
current was measured with a Keithley 6487 amperemeter 
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Figure 1. Experimental setup. The configuration sketciied was chosen for current measurements, being the ampereme- 
ter denoted by -A-. For rate determinations the amperemeter was replaced by two amplifying stages and a linear Fan-in 
Fan-out with one output signal sent to the discriminator and scaler and the other one to the ADC. 






Figure 2. Illustration of typical beam spots in the Polaroid film stack as a function of the number of film. From 
top to bottom: A = 143 ± 14 mm^, A = 36 ± 3.6 mm^, A = 6 ± 0.6 mm^. Left-right shows the increasing direction 
(downstream) of film number. A typical choice for the determination of A is the central column (5). Film 2 is over- 
exposed. Comparison between 5 and 9 shows a spot fading away in a rather uniform manner, suggesting that the 
beam-spot in 5 is highly uniform, as expected. 
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Figure 3. Energy loss spectra obtained witti the s=4 mm anode-pitch MWPC after calibration. The characteristic 
escape peak from Argon is visible under illumination with a Fe^^ source (resolution Oe/E = 15%) but disappears when 
using the X-ray tube. The broader energy distribution in such a case stems from the non-monochromatic nature of the 
bremsstrahlung photons. Single (for the tube) and double (for the source) Gaussian distributions have been fitted to data. 



Once the current and the avalanche rate are measured, the detector gain can be obtained as: 

m=— (2.1) 

where / is the total current, r the avalanche rate and n„ the average number of electrons initially 
released by the X-ray photon. 

For studying the behaviour of the gain as a function of rate, the determination of the later can 
be simplified by resorting to its proportionality with the tube current, that was verified over the 
whole range. Only one direct determination of the rate (obtained at low rates) is therefore needed 
at the beginning of each 'rate scan' and higher values can be obtained by proportionality with the 
tube current. 

In some particular cases Uke for the determination of the characteristic gain curves, where 
the nominal rate was not too high, it was useful to perform direct rate measurements with some 
regularity and care was taken of checking that the electronics dead-time had little impact in the 
results. The dead-time of the electronics system was estimated by resorting to the afore-mentioned 
proportionallity with the tube current: measurements were performed in pure CO2 at m ~ 5 x 10^ 
and a phenomenological model of only one parameter was used to fit the data (Fig. ^, yielding the 
correction formula needed translate the measured rate r^eas to the real rate r. 

r = --^ln{\-AT r„,eas) (2.2) 

with Ar in the range 2-3 /^s, slightly depending on the gain and mixture. Due to this dependency, 
only dead-time corrections not in excess of 20% were allowed in practice. 

A comparison between these two independent methods (proportionality with the tube current 
and dead-time correction) was nevertheless utilized as a cross-check in a number of cases, showing 
good agreement up to moderate rates. In general, the former method (proportionallity with the 
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tube current) was preferred for the measurements of the rate capability and the second (dead-time 
correction) for the determination of the characteristic gain curve. Faster electronics (required for 
the final application in CBM) is currently under development. 




Figure 4. Electronic dead- time as measured in pure CO2 at m ~ 5 x 10^. The points indicate the measured rate as a 
function of the X-ray tube current. The continuous line is a 1 -parameter phenomenological fit introduced in the text. 



All measurements were performed in a monitored (not controlled) atmosphere with the tem- 
perature, pressure and oxygen content showing values in the range T = 22-25 "C, P= 985-1010 
mbar, f^^ = 8-20 ppm. Fresh gas was injected at 2 detector volumes per minute. The voltage 
was corrected to the standard ambient value {To = 20 °C and Po = 1000 mbar) by introducing the 
reduced potential V* = , as is customary. 



3. Gain curve measurements and chamber stability 

The characteristic gain curve is a necessary input when attempting a quantitative description of the 



rate capability of MWPCs [[10|]. On the other hand, the maximum operating voltage before glow 
discharge or breakdown appears in a chamber must be also scrutinized critically. The measurements 
devoted to these aspects are compiled in this section. 



The experimental determination of the gain curve is based on eq. 2.1. The chamber was irra- 
diated with the X-ray tube and, after optimization of the filters and tube voltage, a spectrum with a 
peak energy of Ex-ray — = 6.6 keV and 15 — 30% width (worsening when increasing the fraction of 
quencher) was produced. The main advantage of this procedure, as compared to irradiation with a 
Fe^^ source, is that the very strong decrease in absorption cross-section with the CO2 concentration 
(and therefore in counting rate) can be accommodated by adjusting the intensity of the tube with 
the additional help of 2 Al filters, so that the measurements of the characteristic gain curves can be 
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easily taken at the same avalanche rate. This approach largely increased the reproducibility of the 
measurements under different conditions. 



Once the rate and the current are known, the gain can be calculated from eq. gjj and the initial 
number of electron-ion pairs «o '■ 

(3.1) 

W=[—^ + '-^\ (3.2) 

y Wnohle W^.„J 

being W the effective average energy to produce one electron-ion pair, Wnotie and w^^^ the corre- 
sponding values for pure noble gas and CO2, respectively, and /coj the fraction of CO2. A slight 
shift of the peak energy Ex-ray = 6.6 — 7.1 keV was observed when using a second filter. This was 




taken into account for evaluation of eq. 3.1 



In order not to bring the amperemeter to the limits of its sensitivity, a primary rate r = 120 kHz 
was kept constant for the measurements at low gains (m < 5 x 10'^), yielding a minimum current 
of 5 pA (corresponding to m ~ 1) in pure CO2. The chosen rate represents a compromise between 
a low system dead-time (Fig. ^ and small space-charge, on one hand, and a reasonable current at 
low voltages on the other hand. The uncertaintinty of the current was determined through the rms 
of 30 consecutive measurements, resulting in a statistical dispersion of ~ 1%. At m ~ 1 the current 
stabilized only after several minutes (maximum 30 min), showing an uncertainty of ~ 10%. 

Aiming at reducing the effect of space-charge in the characteristic curves, gains above m ~ 
5x10"* were determined at a reduced rate of r = 10 kHz. After correcting for the system dead-time 
both the high and low gain regimes were found to be in good agreement. No significant difference 
was observed between the gain curve obtained with the Fe^^ source as compared with the tube, for 
a reference Ar-CO2(80-20) mixture. 

For illustration, the gain curves as a function of the reduced anode potential V* are shown 
in Fig. ^ for the case of Xenon mixtures in the = 3 and 4 mm -pitch chambers. The region 
m = 1 is well described when obtaining rio from the tabulated values of Sauli |12], but an excellent 
agreement was also verified for the Ne and Ar mixtures. Overall, the behavior of both chambers 
is very similar, regarding voltage dependence. For chambers of pitch s = 2 mm (not shown) the 
gain curves are substantially different due to the increasing influence of the neighboring wires in 
the field close to the anode. A systematic study of the gain behavior of these chambers will follow 
to this publication. 

Fig. ^ shows an exponential fit in the region m = 10^ - 10^ (solid line), that is used in the next 
section, together with an extrapolation to higher gains (dotted line). The deviation of data from the 
fit was identified as the onset of a self-sustained discharge in the chamber^ . From the crossing point 
of the up-extrapolated gain (dotted line) and a linear down-extrapolation (dot-dashed line) of the 
last two data points one can obtain a crude estimate of the onset of the process. For Xe mixtures 
the chambers were stable up to a maximum gain m,nax = 3 x 10^, being little dependent on the 
fraction of quencher. After prolonged operation in the discharge regime, the chamber required a 
power cycle to restore its previous performances. 



' The chamber continues to draw the same current even when the tube is switched off. After little training, the process 
can be identified by direct observation in the scope. 
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Figure 5. Up: gain m vs reduced voltage V* for the s = 3 mm -pitch chamber operated in Xenon with different fractions 
of CO2. Down: like up-picture but for the s = 4 mm -pitch chamber. The inset zooms in the region where data has been 
fit to an exponential (logjg (m) = 3 — 5) and also the regime of self- sustained discharge (logjQ (ra) > 5.5). 
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4. Rate capability measurements 



4.1 The Mathieson model 



The theoretical description used in this work has been introduced elsewhere [10]. For simplicity, 



an exponential dependence of the gain with the applied voltage has been assumed, such that: 

ln{m)=aV + b (4.1) 

that is well satisfied for our chambers in the region m = 10^ - 10^ and IS a convenient parame- 
terization provided the gain was never reduced below those values, when operated at high photon 
fluxes. Following [p^, the gain behaviour as a function of the flux can be expressed as: 

/ a^{f)no{f) sh?- d,n{s,h) m \ 
V 2jJ.a[f) CL{s,h,ra) \nm-b{f) J 

with qe being the electron charge, the gain at zero rate, no the initial number of electron-ion 
pairs, the mobility of ion-specie a, f the fraction of quencher, Ci the wire capacitance per unit 
length and the avalanche flux (in [T^^L^^]). The correction for finite area of illumination, here 
denoted as dm, is a function defined between and 1, that solely depends on the chamber geometry 
and the shape of the beam spot, under the assumptions of |[To|]^. The functional expression of Cl 
is: 

= w / A (4.3) 



with ra the anode radius and being well approximated in the limit s/h < 2n by: 

—e'"'/' (4.4) 



271 

It becomes apparent that, once the characteristic 'm vs V curve has been measured (a, b are then 
known), the only free parameter for every mixture is Haif)- Therefore, within the Mathieson 
model only one curve 'm vs <p' is needed experimentally for a given mixture (for any nio) and the 



behavior for different nio can be derived from eq. |4-.2| . Moreover, the mobility of the ions is usually 



well described by the Blanc's law ||13|], that in our case simply reads 



J_ = hJ£0^ + Jco^ (4.5) 

Ma l^a, noble l^a,C02 

where IXa,nobie and IXa.cOi the mobilities of the drifting ion a in the corresponding pure gas. 
In case of pure gases the nature of a must change but, if two rate curves are obtained at, say. 



fcoi = 20% and fcoi = 80%, then eq. together with eq. would allow for a complete 
characterization of the rate capability of the chamber for all the possible admixtures of two given 
gases (except maybe in the limit of pure gases). In general different drifting species may coexist. 



but on the basis of a generally accepted fast charge-transfer mechanism [ |12| ] only one drifting ion is 
considered here for each noble gas-C02 mixture, its nature not depending on the gas concentration 
(except maybe in the limit of pure gases). The complex nature of the ion drift and its high sensitivity 



In particular ion diffusion is not included in the model. 
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to the presence of impurities [13] advice, nevertheless, to verify Blanc's law experimentally for 
every particular case. 

The main experimental difficulty involved in the precise evaluation of the rate capability of 
a chamber can be observed after inverting eq. [4.2|: if the rate capability 0^ is defined as the flux 



needed for causing a certain fractional gain drop F , the following expression can be obtained: 

;n((l FK)-^ 2,(/) C.(..,...) ,,,, 

( 1 - F)m„ (/) rio (/) sh?-dm {s,h) 

and making use of the fact that the logarithm is a slow varying function of its argument: 

(t>F F 

Therefore, an experimental uncertainty in the determination of nio by 5% will cause an uncertainty 
in the determination of the fractional drop at F = 10% by 50%, and this uncertainty will be directly 
propagated to the estimated rate capability. 

An illustration of the measured gain behavior as a function of the avalanche flux is shown in 
Fig. I for the case of Xe-CO2(90-10) and Xe-CO2(80-10) in the 5 = 4 mm chamber. The statistical 
uncertainty is smaller than the size of the data points. Indeed, not all the mixtures were equally 
well characterized: in view of the final application of this work (TR detection) the measurements 
in Xenon were performed at 3-4 different values of mo for every concentration of quencher /co,. 
For Argon and Neon -based mixtures almost 2 curves at different nig were measured in average for 
every value of fcoj- The flux has been already extrapolated to the case of uniform irradiation 
as <p dfn<p, taking into account the finiteness of the beam-spot. As shown in the next section, 
the correction factor involved in such an extrapolation dm amounts to ~ 0.5 for the bulk of the 
measurements presented here. 

4.2 Extrapolation to infinite area of illumination 

The 'Mathieson model' includes a prescription on how to correct for finite beam size effects (Hi 



Jllp) and even extraordinary experimental agreement has been found for the 1-D case [ |14| ] (finite 
irradiation along the direction of the wires, infinite in the transverse direction). To the authors' 
knowledge, no dedicated study exists for the 2-D case. 

Following the original formulation of Mathieson, the voltage drop AV{x,y) at the anode plane 
caused by a charge distribution with a certain x-y profile can be obtained from an electrostatic cal- 
culation assuming a constant space-charge density of the drifting ions in the chamber. A correction 
factor can be then introduced as: 

U^,y) = ^ (4.8) 

where AVg is the voltage drop at the anode plane in case of uniform/infinite irradiation. A calcula- 
tion of the factor dm{x,y) for a squared beam profile^ can be accomplished by using expressions in 



[]15[], yielding: 

^'"("'^) = ^X Jo "^^'"^^^ kMkl+kl) 



^The experimental oval-like beam-spot will be assimilated to a square, for the sake of simpler theoretical description. 
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Figure 6. Up: gain as a function of the avalanche flux for the s = 4 mm chamber operated in Xenon with 10% CO2. 
Down: like up-figure but for 20% CO2. The flux has been extrapolated to the case of uniform irradiation. This correction 
implies a scaling-down of the measured flux by a factor 2 for our setup (see text). 



cosh {h{kj + kjy/^) 



(4.9) 



Following the notation of [11], L is the beam size along the wires (x axis) and H across them (y 



axis). When H ^ <x> (1-D limit) dm{x,y) must be averaged in the interval ±L/2 in order to obtain 



the desired correction factor, as has been shown in [ 10 ]. This factor as in eq. ^ would be naturally 
interpreted as an 'effectively' reduced flux d^x <p or, equivalently, as a certain particle rate over an 
'effectively' increased area, as compared to uniform irradiation. Remarkably, according to [ p^ the 
correction factor dm{L,H ,s,h) depends only on the chamber/beam arrangement. 

Despite this apparent simplicity, the 'effective' voltage distribution at the wire positions yic 
(AVe{x,yi^)) seems to be needed in order to generate the exact 2-D distribution AV{x,y) that allows 



for a precise description of data []14|]. A generalization to the 2-D case by averaging dm{x,y) over 
the H /s irradiated wires was introduced in [ [TTI ] based on the magnitude AVe{x,y). The original 
approximation will be used in this work, AV (x, j/t) ~ AVe{x,yic), for the sake of simplicity, so that 



eq. 19 can be used directly. 



In the continuous case where the density of wires is infinite, a direct integration of eq. 
yields the desired average: 

1 rL/2 rH/2 

dm,c = -rn / d,n{x,y)dydx (4.10) 
Ltl J-L/2J-H/2 

while in the discrete case the integral over y must be replaced by a discrete average at the position 
of the wires: 

1 ^ 

^'" = 777/ \dn,{x,yk)dx (4.11) 

NL.J-L/2^^ 

where k runs up to the wires within the area of illumination. If the discretization of space would 
influence the value of d,„ also the relative position chamber-beam will. This was evaluated by 
simulating chamber-beam displacements along y in 0.25mm intervals for different pitches s = 3,4 
mm and a beam-spot L ~ // ~ 6 mm, yielding a deviation from the continuous value dm,c within 
±10%. 

In order to better assess the role of the beam size, a specific set of measurements of the rate 
capability at various beam sizes was accomplished. Taking advantage of the fact that the correc- 
tion factor dni{h,s,L,H) should not depend, within the model, on the gas mixture, its value was 
evaluated from measurements of the rate capability at 10% drop (0io) in pure CO2 (Fig. The 
correction factor was obtained experimentally after normalizing by the rate capability obtained at 
infinite beam size dm = 0io(°°)/0io(^)- The value of 0io(°°) was determined as the average of the 
three rate capabilities measured at the largest beam profiles. 

Fig. ^ shows the behavior of dm both in data (triangles) and model (circles) for a square profile, 
together with a phenomenological 2-parameter fit to a function dm = 1 — exp(— L/Lo') (continuous 
hne). The theoretical value of dm{L) tends asymptotically to 1 although slower than one may 
expect (even for L = 20 mm (» s, h) the required correction factor is still at the level of 20%)^ and, 
conversely, it drops to zero very fast when L ~ as also seen in data. Since data, model and fit 



"^To reduce the correction below 5% would require of a square spot with its side being larger than 10 cm. 
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agree reasonably well for our nominal beam profile L ~ 6mm, a value of = 0.5 was taken for 
extrapolating the measured fluxes to the uniform irradiation case, and a 20% overall uncertainty was 
estimated for the flux calculated following this procedure. This value corresponds approximately 
to the uncertainty of the mean of the three highest points in Fig. ^ that was experimentally defined 
as ^io{°°). 



1.6 




L [mm] 

Figure 7. Comparison between data (triangles) and model (circles) for the correction factor dm = ^io{°°) /<j>io{L). d„ 
has been experimentally defined as the ratio of the flux at 10% gain drop for infinite area of irradiation divided by the 
one measured for a square beam of side L. (/iio(°°) is defined as the average of the three higher points. Data has been 
taken in pure CO2. A 2-parameter fit is shown also as a continuous red line. 



4.3 Fitting procedure 



A simultaneous description of all data based on eqs. ^ and (Mathieson model including 
Blanc's law) is presented in this section. Following the spirit of [ |I3| ] no assumption is done regard- 
ing the nature of the drifting ion that is denoted by a (in Xe), j8 (in Ar) and 5 (in Ne). There are 
therefore 6 free parameters, namely: Aia,C02' Ma.Xe, M/5,C02' M/3.Ar> Ms.coj' ^^S,Ne- A global fit of the 
gain vs flux curves for the 3 mm -pitch chamber (36 data sets) and the 4 mm -pitch chamber (52 
data sets) was therefore attempted by using the 6x2 free parameters introduced above (2 for each 



noble gas-mixture). The transcendental equation 12 was evaluated numerically by interpolation 
after being tabulated first. 

A minimization was performed by assigning equal weight to the data sets for every mixture. 
This means that each residual is normalized by the measured value and each data set is additionally 
weighted by the number of data sets availables per mixture (for instance, the residuals for each data 
set for the Xe-CO2(90-10) mixture in s=4 mm (6 in total) are correspondingly divided by 6, Fig. 
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Figure 8. Up-left: fit of the rate curves for Xe-C02 (90 — 10) obtained for the 5 = 4 mm chamber. Up-right: the same 
as up-left but for 80 — 20 admixture. Center-left: the same as up-left but for s = 3 mm. Center-right: the same as up-left 
but for Ar-C02. Down-left: the same as up-left but for Ne-COa- Down-right: fit of the rate curves for Xe-C02 (10 — 90) 
obtained for the s = 3 mm chamber. Dashed lines represent the theoretical flux at given fractional gain drops from the 
Mathieson formula with parameters obtained from the fit. A total of 88 curves/data sets have been fitted simultaneously. 



^up). This procedure avoids that the fit is dominated by high gains, and compensates for the fact 
that not all gas mixtures are equally well represented in terms of data sets (see app. 
Two main sources of uncertainties were identified: 
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1. One is related to the estimate of nig that is obtained for each data set as the average of the 
3 gains measured at the lowest rates. This procedure works generally very well but in some 
cases there were not enough points for a good determination. Estimates of nig wrong by few 



percent are very critical (eq. 4.7) and, indeed, a typical disperssion of 1% is present in data. 
In order to calculate the error introduced in the fit due to this effect, was globally scaled up 
and down for all data sets in the range ±3% in steps of 1% and the fit performed afterwards. 
The parameters' uncertainties were estimated from the rms of the values obtained from these 
different fits. 

2. A second source of uncertainty comes from the fact that the measurements were taken along 3 
major campaigns and 3 (slightly) different beam sizes were used. Extrapolations to uniform 
irradiation were made according to Fig. ^ for which a systematic error of 20% in the re- 
calculated flux was assigned. 




Figure 9. Compilation of different curves extrapolated to a typical gain nio = 10'*. Figure aims at illustrating the main 
dependences of the rate capability, together with a detailed comparison with the theoretical model (lines). 



The fitted data are compiled in the appendix ^ while some relevant cases are presented in Fig. 
|8[ in order to illustrate the basic dependencies of the rate capability. The curves at fixed fractional 
drop are shown (F = 5%, 10%, 15%, 20%), obtained from evaluation of eq. 4.6 after substitution 
of the parameters of the fit. 

A different attempt to extract dependencies can be performed by comparing at fixed gain. This 



is largely facilitated if the predicted scaling with nio from eq. is used: 



<Pf (mo) _ ln{{l -F)mo)-b 
^,{m'g) \n{{l-F)m'g)-b nio 



(4.12) 
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The second term corrects linearly for the different total charge while the first incorporates the linear 
dependence of the drift velocity with the anode potential. Different dependencies can be therefore 
made more apparent, as shown in Fig. ^. The data sets have been chosen so that a (relatively small) 



extrapolation by using eq. 4. 12 is done for a ratio (m^ — m'^)/mo smaller than 50%. 
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Figure 10. Up-row: rate capability at 10% gain drop for the s = 3mm and s = 4mm pitch chamber at nig = \(f as a 
function of the fraction of noble gas. Data correspond to Xenon (triangles), Argon (squares) and Neon (circles) based 
mixtures. Open symbols indicate that the plotted value has been extrapolated over more than 50% difference with respect 
to mo- Lines show the Mathieson model with the mobilities for the drifting ions taken from a global fit to 88 data sets. 
Low-row: as up-row but for mg = 5x 10^. 



A more systematic approach can be still devised: the flux for a 10% fractional gain drop from 
a nominal gain nio ('(pio at m,/) can be experimentally determined by extrapolation from the closest 
value measured at m[„ according to eq. i.l2\ The main uncertainty in the determination of 0io is 



the uncertainty in the measurement of rrio, that has been assumed to be 1%. By resorting to eq. 4.7 
an uncertainty of 10% in has been estimated. 

The resulting points can be plotted as a function of the fraction of noble gas, as shown in 
Fig. |l^ for Xenon (triangles). Argon (squares) and Neon (circles) mixtures, for different pitches 
(s) and initial gains (mo). Full points indicate that the extrapolation has been done from a curve 
measured with a real gain m'^ within less than 50% deviation from the extrapolated m„. The lines 
represent the fit to the 'Mathieson model with Blanc's formula'. For convenience, the dependence 
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of the gain parameters a and b with /coj has been parameterized by using 3*^^ order polynomials. 

The fitted mobilities used as input for the model are compiled in table [IJ with their statistical 
and systematic uncertainties. The later (20%) arise from the uncertainty in the flux determination 



due to the beam size correction, that is reflected Unearly in the measured mobilities (eq. [4.6| ). This 
uncertainty is a correlated one for all the measurements. The following features can be observed: 
i) the values of fitted mobilities in the s = 3 and 4 mm chambers are compatible within lc7 with 
exception of HajCe and jU/j co^ ii) the mobilities in noble gases differ in a ratio ~ 1 : 3 : 4 for 
Xe:Ar:Ne. iii) When the results from the two chambers are combined, ion mobilities in Xe, Ar and 
Ne differ with a statistical significance of more than 2a while all the mobilities in pure CO2 are 
statistically compatible within la, iv) the values of IJ.a,Xe are sUghtly higher than expected for Xe+ 
but slightly smaller than other indirect measurements as [0], v) the mobilities of ion j8 in Ar and 
CO2 agree well within errors with the COj hypothesis, as has been measured for similar chambers 



[ 13], vi) ion 5 shows a mobility much closer to the value of Ne+ than CO2 (factor two difference), 
vii) the average is 0.99, providing statistical support to the approach followed in this analysis. 



A compilation of these mobilities can be found in []16[]. 



mobility [mm^V ^s ^] 


s = 3 mm (fit) 


s = 4 mm (fit) 


mean 






155 ±37 ±31 


124 ±16 ±25 


129 ±15 ±26 


0.59 


IJ-a,Xe 


70±6±14 


95 ±13 ±19 


74±5±15 


3 




173 ±18 ±35 


128 ±24 ±25 


157±14±31 


2.25 




214 ±20 ±43 


207 ±30 ±41 


212±17±42 


0.0377 


I^S.C02 


126 ±23 ±25 


129 ±25 ±26 


127 ±17 ±25 


0.0081 


IJ-S,Ne 


298 ±21 ±60 


289 ±30 ±58 


295 ±17 ±59 


0.0604 



Table 1. Ion mobilities from the global fit of the data to the Mathieson model. The first errors are statistical, the second 
ones systematic. 



4.4 Calculations for minimum ionizing particles 

In the CBM experiment at FAIR, most of the particles over the anticipated MWPC-based TRD will 
be minimum-ionizing particles (mips) [^]. It is therefore important to estimate the rate capability 
in such an environment. A systematic survey of energy loss in Xe-C02 has been recently carried 



out over a broad dynamic range [18] where a detailed comparison between data and simulation 
was presented. A calculation purely based on the Bethe-Bloch formula overestimates the measured 
energy loss by some 30%. Hence, for simplicity, the Bethe-Bloch prescription as recommended 
for mixtures [ p^ was used here, but with a normalization taken to reproduce the value dE/dx = 5 



keV/cm for p = 1 GeV pions in Xe-C02(85-15) reported in Q18|]. When going from pure Xenon to 
pure CO2 such a procedure predicts a reduction in energy loss by a factor of 2, due to the smaller 
density. On the other hand, in the above-mentioned reference gas mixture the energy lost by mips 
is approximately 3 keV in the total gas thickness 2h = 6 mm, to be compared with our X-ray energy 

Ex-ray = 6.6 kcV. 

The compilation of the measured rate capability in Xenon mixtures for the s = 3 and 4 mm 
pitch chamber, extrapolated to mips by using the Bethe-Bloch formula obtained as explained above. 
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Figure 11. Up: rate capability at 10% gain drop for the s = 3 mm (diamonds) and s = 4 mm (stars) chambers 
extrapolated to minimum ionizing particles. The lines show the Mathieson model without fit. The extrapolation to the 
.5 = 2 mm case is also indicated (dot-dashed line), by using the mobilities of the i = 3 mm chamber from table in text 
(the necessary parameters of the characteristic gain curve are taken from measurements). The dotted curves (scale on 
the right) correspond to the absorption probability (Pahs) of X-rays of £ = 6 keV. Open points indicate that the data has 
been extrapolated to ra,, from a measured value m'„ different by more than 50%. 



together with the model and an extrapolation to s = 2 mm is shown in Fig. ^ for a typical operating 
gain of nio = 10'* (extrapolation to s = 2 mm indeed requires to use the gain curves measured in 
such a case, that will be published elsewhere). 

The improvement of the rate capability with decreasing s has two main origins: first, the rate 
per wire is proportionally smaller in such a case; second, the higher wire density implies that higher 
fields are required in order to get similar gains, meaning that the coefficient a becomes smaller (eq. 
4. 1 ). The fact that the gain curve is 'softer' in the sense of a being smaller makes the chamber more 
resistant to the space-charge of the ions {^f ~ 1 /a^). The capacitance per unit length Cl has little 
impact for normal values of s. 

On the other hand, the behavior of the rate capability with increasing fraction of CO2 has four 
main dependencies: i) The energy required to create a pair in CO2 is higher than in noble gases 
(exception is Neon, where both are comparable). Still, this fact can account for a 30% increase 
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in Xe-C02 mixtures due to the lower initial number of electron-ion pairs when increasing the 
concentration of CO2. ii) CO2 is less dense than the mixture by a factor amounting by maximum 
2 in Xenon-based mixtures, increasing the rate capability accordingly when going from pure Xe 
to pure CO2. iii) The curves become 'softer' when the fraction of CO2 is increased, resulting in 
an increase by a factor 2-3 from pure Xe to pure CO2. iv) The mobility of the ions can vary also 
sizably between extreme cases: in case of Ne this implies a reduction of up to a factor 3 when 
increasing fcoi while in Xe it increases by a modest 50% when going from pure Xe to pure CO2. 
These effects compete in some cases, as in Neon, where the increase in rate capability due to the 
much increased mobilities at high Neon concentration is compensated by a much harder gain curve. 
This yields a largely flat dependence of the rate capability for X-ray photons as a function of CO2 
concentration (Fig.[l0|). In Xe mixtures all effects i-iv go in the same direction, and up to one order 
of magnitude can be gained in rate capability for minimum ionizing particles, between pure CO2 
and pure Xe, as shown in Fig. [TTI-up. 

The present model has been used to describe recent mips data from M. Petris et al. [p^], 
showing a good agreement. Nevertheless, we believe that the very small deterioration of the per- 
formances observed by the authors does not add extra support to the message conveyed in the 
present work. 



5. Discussion 



Considerable attention has been devoted in this work to the 'beam-size' correction. Nevertheless, 
if the model assumption that 'the correction factor dm depends only on the chamber/beam arrange- 
ment' would hold true one could eventually resort, in general, to the weU measured CO2 mobility 
in CO2 to estimate the factor <i„, by imposing agreement with the Mathieson formula 4.2 in the 
limit of pure CO2 (see [20]). Instead, if one would resort to the theoretical model to calculate 
it, a correction factor that would include i) an accurate knowledge of the beam intensity profile 
and ii) the relative beam/chamber alignment would be needed in order to obtain a precise (yet 
model-dependent) correction. Using the known COj mobility in CO2 as a reference or, alterna- 
tively, taking one single measurement with defocused beam should be probably more sound than 
attempting a detailed calculation. Overall, we consider the results of table |l| as a strong indication 
that reliable predictions can be done for the rate capability of MWPCs once the proper geometric 
corrections are performed and the nature of the drifting ion is known (and its mobility has been 
measured before). Within the precision of the present approach, Blanc's law works well. There 
is, nevertheless, a tendency of the Argon and Xenon mixtures to behave similarly at high gains 
{nio = 5 X 10^), a fact that cannot be accomodated in the model. 

With all the above information at hand, we are in position for discussing the more conve- 
nient MWPC geometry for the envisaged application in CBM. Due to the nature of TR, hadron 
identification/suppression worsens at low momenta where the particle deviates from the minimum 
ionizing regime. Additionally, the nature of the CBM spectrometer makes the detection of high 
energy electrons specially relevant in the context of and detection, while the low momenta 
regime is well covered by the anticipated RICH detector and the TOP wall. Based on that we 
take as a working number the average TR energy released by a = 2 GeV electron in a typical 
CBM-radiator, that is approximately f'j,^ ~ 6 keV [Bl]. This roughly corresponds to the minimum 



-18- 



momentum for which the CBM TR system will provide sufficient ejn suppression, and that we 
consider as our benchmark here. From the tabulated X-ray absorption in Xenon a mean free path 
of Xy(E = 6keV) = 2.7 mm can be readily obtained, being the probability to absorb the photon 
in the gas given by Pats = 1 — exp(— 2/j/A). High absorption probability must be balanced with a 
moderate rate capability, that has a steep dependence with the gap size (eq. 4.6). 



Xe/CO^iSS/IS), s=h 




h[mm] 



Figure 12. Extrapolated rate capability of a MWPC operated in a Xe-C02(85-15) mixture as a function of the chamber 
gap h for mips (axis on left) when the average signal amplitude in number of electrons is fixed to 10^ e^ and i- = h. The 
absorption probability of TR photons for the CBM benchmark case is shown with axis on the right (dashed line). The 
shadowed area represents a comfortable situation in terms of rate capability and photo-absorption probability. Mobilities 
have been taken from the weighted average of the mobilities in the s = 3 mm and s = 4 mm chambers, as obtained from 
the global fit to the Mathieson model. 



A meaningful determination of the rate capability as a function of h can be accomplished by 
making some further assumptions: first, in order to avoid potentially harming operating voltages, 
the ratio s/h is kept constant. Second, the operating conditions are assumed to be such that the 
average number of electrons after multiplication is constant and equal to rie ~ lO^e^ for mips 
(changes in h lead to changes in the initial ionization, that must/can be compensated by re-adjusting 
the gain, so that the induced signal is roughly the same). It is also possible to directly obtain the 
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rate capability dependence from a direct evaluation of eq. 16, when recalling the proportionallity 
of rig with h, but at the price of a reduced signal for small h. In either one case or the other the 
implicit assumption that the gain curve is not strongly depending on h must be made, and only then 
an actual l/h^ behavior (slightly modified by Cl) can be obtained for the rate capability. 

Fig. [I^ shows the expected behavior of the rate capability for MWPCs operated in Xe-C02(85- 
15) as a function of the anode-cathode gap h under the former assumptions, together with the 
absorption probability for our benchmark case in CBM (Ej.^ ~ 6 keV). Indeed, a very narrow set of 
values can be considered as satisfactory (h = [2.5 — 5.5 mm]), with the minimum value determined 
by the condition Pabs > 80% and the maximum by > 100 kHz/cm. 

Being the gap dependence so critical for MWPCs in terms of photo-absorption probability, it 
is worth discussing the possibility of using a drift region (of length D). For typical pad sizes of 1 
cm^ [^, the maximum CBM fluxes would imply a particle load of 1/10 jxs per channel causing a 
10% pile-up probability within 1 /xs. Taking the measured drift times in even in a small drift 
region corresponding to D = 2h, the electron drift would extend up to a typical time of that order 
[Q], making accurate charge measurements very difficult. 

On the basis of the above arguments, staying at a MWPC configuration seems to be the most 
convenient solution for the CBM TR-detector. The possibility of using a mirrored configuration 
with 3 cathode planes deserve further consideration, since it increases the X-ray detection probabil- 
ity while keeping the rate capability. If this improvement is worth the extra mechanical complexity 
or not, will be clarified when the first real size prototypes are built. 



6. Conclusions 

We conducted systematic measurements of thin MWPCs (h = 3 mm, s = 3,4 mm) filled with 
Xenon-C02, Argon-COa and Neon-C02 mixtures under different fractions of quencher and at dif- 
ferent gains. Employing the Mathieson model (including finite beam corrections) and the additive 
Blanc's law for ion mobilities the measured rate capability can be well described theoretically. 

When aiming at efficient TR detection at high rates there is a narrow parameter-space where 
MWPCs are advantageous. The reason for that is the steep dependence of the rate capability 
with the chamber gap /j as ~ 1 //j^ opposed to the exponential behavior of the X-ray absorption 
probability Pats = 1 — exp(— 2/j/A). We found a good compromise in the value h = 3 mm where we 
can comfortably hold rates above 100 kHz/cm^ as demanded by the TRD of the CBM experiment, 
while keeping a typical X-ray absorption probability above 80%. 
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A. Results of the fit 



A compilation of data for 36 mixtures (from a total of 44), fitted according to the procedure in- 
troduced in section 43, is presented here. The 8 missing plots are not shown (but included in the 
fit) for easier graphical representation. Data for each pitch and noble gas mixture (for all measured 
concentrations of the two gases) have been fit with 2 free parameters, namely the mobilities of the 
drifting ion in the pure gases. Due to the main scope of this work, data from Xenon at different 
'low rate' gains are more abundant than from Argon and, specially. Neon. For the later, only 
1-2 curves are taken per mixture, a choice that relies on the trivial scaling of m„ predicted by the 
Mathieson model (eq. 4.12). Each mixture enters in the fit with the same weight, meaning that the 
residuals for each curve are divided by the number of curves per mixture (for instance, the residuals 
in Xe-CO2(90-10) are divided by 6). The residuals are additionally divided by the measured value, 
in order that all data sets enter with the same weight for a particular mixture. 
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Figure 13. Up-left: 2-parameter fit of the rate curves for Xe-C02 mixtures obtained with the s = 4 mm chamber. 
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Figure 14. Up-left: 2-parameter fit of the rate curves for Xe-C02 mixtures obtained with the s = 3 mm chamber. 



-23- 




Figure 15. Up-left: 2-parameter fit of the rate curves for Ar-C02 mixtures obtained with the s = 4 mm chamber. 
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Figure 16. Up-left: 2-parameter fit of the rate curves for Ar-C02 mixtures obtained with the i = 3 mm chamber. 
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Figure 17. Up-left: 2-parameter fit of the rate curves for Ne-C02 mixtures obtained with the s = 4 mm chamber. 
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Figure 18. Up-left: 2-parameter fit of the rate curves for Ne-C02 mixtures obtained witii the ^ = 3 mm chamber. 
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